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Abstract. We carry on a new analysis of the sample of MACHO microlensing candidates towards the LMC. Our 
main purpose is to determine the lens population to which the events may belong. We give particular emphasis to 
the possibility of characterizing the Milky Way dark matter halo population with respect to the LMC one. Indeed, 
we show that only a fraction of the events have characteristics that match those expected for lenses belonging 
to the MACHO population of the Milky Way halo. This suggests that this component cannot explain all the 
candidates. Accordingly, we challenge the view that the dark matter halo fraction of both the Milky Way and 
the LMC halos are equal, and indeed we show that, for a MACHO mass in the range 0.1-0.3 Mq, the LMC halo 
fraction can be significantly larger than the Milky Way one. In this perspective, our main conclusion is that up 
to about half of the observed events could be attributed to the LMC MACHO dark matter halo. 
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1. Introduction 



Since the original proposal of lPaczvhskil l|l986f) . microlens- 
ing has shown to be a powerful tool for the investigation 
of the MACHO composition of the galactic halos. The mi- 
crolensing surveys towards the LMC and the SMC probed 
the existence of compact halo objects along these lines of 
sight, however, the assessment of these results with re- 
spect to the actual fraction of the Galactic halo in form 
of MACHOs is still highly debated. 

The MACHO collaboration reported the detection 
of 13-17 microlens ing candidates towards the LMC 
llAlcock et al . 2000b), arguing in favour of a MACHO halo 
fraction of ~ 20% of objects of ~ 0.5 Mq, and estimating 
a microlensing optical depth towards the LMC of r = 
1.2+Q3 x 10 -7 . The reported microlensing rate towards 
the LMC significantly exceeds the expected one from 
known visible components of our Ga laxy. Further analysi s 
mainly confirmed these conclusions l|Bennett et alJl2005|h 
On the other hand, the EROS collaboration, out of obser- 
vations towards both the LMC and the SMC, put more 
and more lower upper limits on the MACHO contribution 
to galactic halos jLasser re et al.ll2000t lAfonso et al.l [2003: 
iTisserand fc Milsztainll2005|h that are no longer compat- 
ible with the MACHO results. 

These questions have been addressed also by mi- 
crolensing surveys towards M31 l|Ansari et al.lll997|) . and 



first evidences for a MACHO contribution along this lin e 
of sight have been re ported llCalchi Novati et al.l 12005^ . 
although challenged in lde JongeTlibTl)2006|) . Overall, the 
picture remains unclear, in particular the problem of the 
nature and the location of the observed events. For ex- 
ample, [Greer^^Jedamzik| l|2002() have cast doubts on the 
interpretation of the microlensing data towards LMC as 
due to a dominant lens population made by MW halo 
MACHOs. Indeed they show, at high level of confidence, 
that the distribution of the duration of the observed mi- 
crolensing events is significantly narrower than what is 
expected from a standard halo lens population. 

We have addres sed s ome of these questions already in 
IJetzer et all l|2002|) and iMancini etafl l|2004j) (hereafter 
respectively Paper I and Paper II). A main issue in both 
works is that the microlensing events towards the LMC, 
observed by the MACHO collaboration, do not belong nec- 
essarily all to the same lens population. In particular in 
Paper II we have co nsidered the issue of self le nsing in 
the framework of the Ivan der Marel et al.l l|2002j) picture 
of the luminous components of the LMC. The main con- 
clusion, based both on the predicted number and charac- 
teristics of self-lensing events, was that self lensing alone 
cannot explain all the observed events. 

In the present work we extend the analysis by fully 
considering the LMC and MW dark matter halo MACHO 
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lens populations. The main question we want to address is 
whether and to which extent events due to the lens popu- 
lation residing in the LMC halo can contribute to the ob- 
se rved rates. T his issue was first proposed and discussed 
in iGouldl (1993). Here we consider again the question tak- 
ing into account the results of the MACHO collaboration 
and the most recent results on the modeling of the LMC. 

The paper is organised as follows. In Sect. Owe review 
the models of the LMC and MW components we consider. 
In Sect. |31 we discuss the evaluation of the microlensing 
rate and present the results for the expected number and 
duration of the microlensing events. In Sect. 01 we carry 
out our new analysis and present our main conclusions on 
the LMC MACHO halo contribution, and in Sect. we 
present our conclusions. 

2. Models 

We consider the LMC as composed by a luminous part, the 
bar and the disc, plus a stellar and a dark matter halo. We 
also include the MW dark matter halo, but we do not in- 
clude the MW disc and spheroid populations. These com- 
ponents have already been shown t o give smaller contr i- 
butions than the LMC self lensing l)Alcock et al.ll2000b|L 
which we instead include. Accordingly, we exclude from 
our analysis the only one event whose lens, upon a direct 
search, has been acknowledged to be part of these compo- 
nents. 

For the structure and dynamics of the luminous com- 
ponents of the LMC we follow closely P aper II. Out of the 
analysis of Ivan der Marel et al who derive their 

results on the assumption that the carbon star population 
is representative of the bulk of the LMC disc stars, we 
take up the results on the LMC disc intrinsic ellipticity, 
vertical thickness, line-of-sight velocity dispersion and ro- 
tation curve and the coincidence of the centre of mass of 
the disc and the bar components. We assume the following 
values for the bar and disc components Mw + M disc = 
(2.7 ± 0.6) x 10 9 M^ llvan der Marel et al.112002^ . and 
-Wbar = 1/5-Mdisc l|Gvuk et al.ll2000(K that we consider to 
be both centered at a, S = 5 h 27.6 m ±3.9 m , -69.87°±0.41° 
J2000) at a distance fro m us of Dq = 50.1 ± 2.5 kpc 
van der Marel et alJl2002|) . We use the same density star 
distribution than in Paper II, characterized by a vertical 
distribution for the exponential disc described by a sech 2 
function and a boxy bar, with a gaussian profile along 
the major axis and the section described by a exp(— r 4 ) 
function. The vertical structure of the LMC has been re- 
cently the object o f great debate (for a discussion see 
van der Mare]H2o"f)l . We have considered this issue, with 
respect to the expected self lensing signal, in Paper II. In 
the present analysis, where we focus on the contribution 
of the two halos, we do not enter such discussion and con- 
sider only the configuration with coplanar disc and bar. 

The presence of a significa nt LMC stellar halo pop- 
ulation is a matter of de bate ijMinniti et al.ll2003l lAlvesl 
l2004t iGallart et al.ll20oj) . In the present analysis we in- 
clude the contribution of such a possible component fol- 



lowing lAlvesl ( 2004) who proposes a spherically symmetric 
spheroid with density profile 

P = Po{l + ^j , (1) 

with k = 3/2, central density po = 6.3 x lO 6 M0kpc -3 
and core radius etc = 1-42 kpc for a total mass, within 8.9 
kpc, of 0.35 x 10 9 M Q somewhat smaller than that of the 
bar. The optical depth profile of this component shows 
a near-far asymmetry due to the LMC disc inclination 
whose overall shape recalls that of the optical depth profile 
of the LMC MACHO halo component (Paper II, Figure 
4), with a maximum value around 0.9 x 10 -8 reached in 
corresponden ce of the field MACHO 6. 

Following lvan der Marel et alJ l|2002() , who present ob- 
servational evidences based on the rotation curves, we in- 
clude a significant LMC dark matter halo component. We 
assume a total LMC mass of 8.7 x 1 9 Mfr> within 8.9 kpc 
with a truncation radius of 15 kpc ijvan der Marel et alJ 
2002), that we consider to include also the contribution 
from the stellar halo. We assume a spherical isothermal 
model (k = 1 in Eg. ^) with a core radius of etc = 2 kpc 
l|Alcock et al.ll2000bl) . 

This spherical symmetric configuration might not be a 
realistic description of this component given that the dy- 
namical environment of the LMC can induce tidal distor- 
tions and disruptions especially in the outer parts. To take 
into account this issue, in Paper II we have compared the 
LMC halo optical depth profiles for a spherical and an el- 
liptical configuration (Fig. 4 and 5 respectively): the over- 
all shape is similar even if in the latter case the maximum 
value rises by about 20% and the near-far asymmetry is 
enhanced. However, as we lack any strong constraint, we 
prefer not to introduce a further parameter in the present 
analysis, therefore we consider only the spherical configu- 
ration, in view also of the possibility to carry out a more 
direct comparison with previous works. 

In this same perspective for the MW dark matter halo 
we consider the "standard" isothermal profile with a core 
radius of 5 kpc, local density 7.9 x 10 6 M kp_c_^_anda 
distan ce from the Galactic centre of 8.5 kpc l|Alcock et alJ 
l2000bft . 

3. The microlensing rate 

The main tool of investigation we use is jy-, the dif- 
ferential rate of mi crolensing events with r espect to the 
Einstein time TV. jde Ruiula et all Il99li iGriesd Il99li 
iRoulet fc Molleracbl fl997). This allows us to make predic- 
tions on the timescale, the number and the spatial distri- 
bution of the expected events, that we can compare with 
the corresponding observed quantities. With respect to 
the self-lensing configuration we have analysed in Paper 
II (section 4.2), we can no longer adopt the useful approx- 
imation = x rs 1. Moreover we now have to take fully 
into account the bulk velocity of the LMC components 
and the relative motion b etween the LMC and the MW 
I van der Marel et aljEoO^ . 
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The source stars belong to the luminous components, 
disc or bar, of the LMC, whereas the lenses can belong ei- 
ther to the LMC or the MW halo. We assume an isotropic 
Maxwellian profile 1 for the velocity distribution for both 
lenses and sources. For the flattened LMC luminous com- 
ponents this is a rough approximation, still, it gives a fair 
description of the average properties of these populations, 
that we consider to be sufficient in the present framework. 

For the LMC disc component we cons i der th e ro- 
tational velocity as in Ivan der Marel et alJ l|2002h with 
a = 20.2 km/s. For the LMC bar stars we use a larger 
value of th e velocity disper sion than for the disc, a = 
24.7km/s (|Cole et alJ l2005|) . For both the LMC halo 
components, stellar and MACHO, we use a — 46 km/s 
l|van der Marel et alJl20o3 lAlvesl l2004|) (we have tested 
that our results remain qualitatively unaltered by chang- 
ing this value up to 20%); for the MW halo, a — 155 km/s. 

The expression for the random motion velocity for the 
lenses reads 2 



■ XV R 



(2) 



where vi± is the velocity relative to the microlensing tube 
at position x, whose modulus is the ratio between the 
Einstein radius and the Einstein time (i>i± = Re/Te), 
v s j_ the random component of the velocity of the sources. 
All the bulk motions are included in A i , defined as 



A±_ = V Q X + X («LMC1 - *OJ 
V («LMCi + l>2,drift_L) 



f s , drift ±) 



(3) 



where r\ — 0, 1 for lenses in the Galaxy and in the LMC 
respectively, vz,drift_i_ is the drift velocity of the lens star 
belonging to the LMC disc (for the bar sources as well 
as for the halo lenses we only consider a random motion 
component) and a tilde over a vector indicates a quantity 
measured by an observer comoving with the MW centre. 
In the self-lensing configuration, A± 0, and Eq. 
reduces to Eq. (13) of Paper II. 

We call a the angle between the inner normal to the 
tube, n, and the source velocity v s ±; 9 the angle between 
h and Vi±, 9 6 (— n/2, n/2); ip the angle between v s ± and 
A±; so that the angle between vi± and v s ± is a — 9, and 
that between vi± and Aj_ is 7 = a + p — 9. Both a and 
(p vary in the range (0, 2tt) (Fig.^l. 

In the case of A± = all the integrations over the an- 
gular variables can be carried out analytically. This is now 
only possible for the 9 variable. Altogether, after an ana- 
lytical integration on the modulus of the source velocity, 
we are left with the following expression of the differential 



1 The Maxwellian profile of the velocity distribution is the 
first term of a series expansion in terms of Gauss-He rmite 
moments Jvan der Marel fc Framll Il993t iGerhardl Il993f) . See 
Section 3.2 of Paper I. 

2 The velocity components parallel to the line of sight are in- 
tegrated out, and the subscript _L indicates the vectorial com- 
ponent in the plane orthogonal to the line of sight. 




Fig. 1. Scheme of a section of the microlensing tube with 
indicated the positions of the vectors and the angles in- 
volved. 



rate with respect to the Einstein time Te'- 
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x 2 of (vi± sin(a) — As'm(ip)) 2 



1 - Erf 



a 2 2(„ 2 + x 2 a 2 ) 
xa s (vi± cos(a) 



A cos(p)) 



(4) 



where the normalization factor M is the integral over the 
line of sight of the sources 



dD os p s (D os ) , 



having assumed that the number of detectable stars varies 
with the distance as D~ 2 . We define x m - m = d min /D os , 
A = \A\; the integration limits along the line of sight, d m i n 
and G? max , represent the distances from the observer to the 
intersection with the LMC tidal surface, for lenses in the 
LMC, whereas d m i n = for lenses in the Galactic halo. We 
use solar mass units, defining /i = j^-, where M is the lens 
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mass. For lenses belonging either to the disc or the bar of 
the LMC, as well as for the stellar LM C halo contribut ion, 
we use the exponential function as in IChabrierl fcOOlft for 
the mass function with integration limits /i m i n = 0.08 
and ^ max = 1.5. If the lenses belong to the MW halo or 
to the LMC halo, the would be MACHOs, we use a set of 
delta function with values going from 10 -5 up to 10 M Q . 

Finally, to compare with the results of the observa- 
tions we have to take into account the expression for the 
detection efficiency, so that we obtain 



dT \ ~ ^ FIT ) 

did -d^' f(TE) - 



(5) 



To discuss the results of the MACHO col laboration we 
use the efficiency as a function of Te as in lAlcock et alJ 
2000b) for which an analytical approximation is given in 



Paper II. We take in to account the correction reported in 
iBennett et all <|2005h . 



3.1. Expected number and duration 

Following the prescription outlined in the previous section, 
we evaluate the differential microlensing rate, for each dif- 
ferent lens population we consider, along the lines of sight 
towards the MACHO fields. In particular for each lens 
population we calculate the expected number of events 
per field as 



(6) 



r°° dr 

Afield = -Eficid / -j^r £{T E )dT E , 
Jo dT E 

where t he "field exposure" -Efieid is defined in lAlcock et alJ 
d2000bh . as the product of the number of distinct light 
curve per field and the relevant time span. Furthermore, 
we pay attention to eliminate the field overlaps. 

To characterize the expected timescale we report the 
median value of the asymmetric distribution (dT/dT E ) e , 
together with the values ?e 16 % and T Ey84 % that single 
out the extremes of the 68% probability range around the 
median. 

In TableHlwe report, for both MW and LMC MACHO 
lens populations, the total number of the expected events 
in all the MACHO fields, in the case of a full MACHO 
halo, together with the expected timescales. We recall that 
the expected timescale varies with the square root of the 
mass of the MACHOs, and that very short timescales are 
strongly suppressed by the detection efficiency function. 

The predicted durations turn out to be almost inde- 
pendent from the position for both the halo populations 
we consider, whereas the issue of the variation of the ex- 
pected timescales with the position across the fields for 
the self-lensing population has been discussed thoroughly 
in Paper II. Indeed, for lenses in the MW halo we find a 
dispersion of the median timescales towards the different 
fields smaller than 1%. In the case of the LMC lenses the 
dispersion is only slightly larger, at most ~ 5%. 

The expected total number of events due to the LMC 
stellar halo turns out to be ~ 0.6, about half of those due 



lens mass 


MW 




LMC 




M 


T E (days) 


iVoxp 


Te (days) 


Nexp 


1CT 5 


o n+3.U 
J *°— 1.3 


0.9 


o n+a.U 
J -°-1.3 


0.4 


1CT 4 




8.3 


3.5l?; 3 


3.2 


itr 3 




52 




13.5 


icr 2 




115 


1311! 


17.3 


0.1 


20+ 15 


97 


361- 


9.9 


0.2 


ZD -10 


82 


48+^ 


7.6 


0.5 


4 1 +29 
^-16 


59 


? r+46 
'°-27 


5.0 


1 


55t 3 2 ? 


44 


103+H 


3.4 


10 


1491™ 


12 


245lg* 


0.5 



Table 1. Expected duration, median values with 68% 
CL errors, and expected number of events for a full dark 
matter halo, respectively averaged and summed over the 
MACHO fields, as a function of the MACHO mass. 



to the LMC disc-bar self-lensing contribution (Paper II). 
The expected median timescale, averaged on the 30 fields, 

+43 



is T E = 45I23 days- 
Overall, we recover the result l)Alcock et al.1 EoOOb) 
that stellar lensing alone cannot explain the signal, so that 
most of the detected events must belong either to the MW 
or to the LMC dark matter halo. 

4. The LMC MACHO contribution to 
microlensing events 

A straightforward conclusion to be drawn out of the re- 
sults on the expected number of events due to the dark 
matter MW and LMC halos is that, with the implicit hy- 
pothesis that the halo fractions in both the MW and the 
LMC halos are the same, most of the lenses should indeed 
belong to the MW halo. Our aim is to challenge this point 
of view. 

First, we recall the current status about the microlens- 
ing events observed by the MACHO collaboration. Next, 
we carry out a statistical analysis of the observed charac- 
teristics of the events (timescale and spatial distribution) . 
The purpose is to determine to which extent the available 
data allow to distinguish between the two halo popula- 
tions. Eventually, using a likelihood analysis based on the 
microlensing rate, we study whether a viable solution to 
the MACHO puzzle can come from a significant contribu- 
tion of a lens population belonging to the LMC halo. 

4.1. The microlensing MACHO candidates 

In the final analysis of 5.7 years of data in 30 fields to- 
wards the LMC, the MACHO group presented two sets of 
microlensing candidates, set A and set B, with 13 and 17 
candidates respectively, the for mer being a subsample of 
the latter l|Alcock et alJl2000b|) . 

Further works allowed to get more information on 
some of these candidates. The len s for the event LMC - 
5 is located in the Galactic disc l Alcock et al.1 EoOlrJl . 
LMC-22 has been identified to be very likely a supernova 
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l|Alcock et al.ll2001a|) - LMC- 23 has been acknowl edged to 
be a probable variable stars {Bennett et al.ll2005l). LMC-9 
is a do uble lens system with caustic crossing (jAlcock et alJ 
l2000a|) . The microlensing candidates LMC-9, LMC-20, 
LMC-22 and LMC-27 belong to the set B only. For most 
of the remaining events a further photometric follow-up 
allowe d to confirm the mic rolensing origin of the flux vari- 
ation {Bennett et alJl2005|) . 

In the present analysis we restrict to an homogeneous 
set of Paczyhski-like events, therefore we exclude LMC-9. 
Furthermore, we do not include the Galactic disc lens pop- 
ulation, so that we exclude LMC-5, as well as all those can- 
didates whose microlensing origin has been put in doubt. 
Accordingly, in the following we shall consider a subset 
of 13 events taken from the original larger set B, from 
which we exclude the candidates LMC-5, LMC-9, LMC- 
22 and LMC-23. Furthermore, we have verified that our 
main conclusions would not change had we started from 
the smaller set of 11 events subsample of the original set 
A, just excluding the candidates LMC-5 and LMC-23. 

4.2. Duration and position: a statistical analysis 

As previously noted, the expected timescale distributions 
for microlensing events due to lenses either in the LMC or 
in the MW halo are almost independent from the position. 
This property allows us to carry on an analysis in which we 
compare the observed timescales to the predicted one for 
each population. In particular we investigate whether it 
is possible to draw from such an analysis any conclusion 
on the relative fraction of the MW over the LMC dark 
matter halo events. Here we neglect the stellar lensing 
contributions. 
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Table 2. Kolmogorov-Smirnov test results 



LMC and the MW halo, in correspondence of two values of 
the MACHO mass (0.2 and 0.5 M Q ). Superimposed, the 
vertical lines indicate the Einstein time of the observed 
events. To further investigate this issue, we make use of 
the Kolmogorov-Smirnov test (hereafter KS). This allows 
us to test the null hypothesis that the events are drawn 
from a given population. The resulting KS coefficient gives 
the significance level of the test. In the first place, we apply 
the KS hypothesis test separately to the two populations 
of lenses in the MW and the LMC halo. Then we intro- 
duce a parameter a, defined as the ratio of the MW events 
over the total (MW plus LMC events), in order to explore 
the possibility that an intermediate solution, with the two 
populations mixed, has to be preferred. To this purpose we 
look for the value of a that minimises the corresponding 
KS coefficient. 




100 150 
IE (days) 



Fig. 2. Normalized differential rate (dT/dTj^) e for both 
MW and LMC halos for 0.2 and 0.5 M , dashed and solid 
lines respectively . Superimposed the value of the observed 
durations. The y-axis values are in 10~ 2 units. 



For a given value of the MACHO mass, the expected 
LMC median timescales are larger than the MW ones 
l|CTOiildlll99,i . In Figure we show in the same plot the 
normalized, differential rate distribution for lenses in the 
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Fig. 3. Kolmogorov-Smirnov and a coefficients as a func- 
tion of the MACHO mass, a is the ratio of the MW over 
the total number of dark matter halo events. 
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To take into account the variations of the microlensing 
rate across the fields, mainly in the case of lenses in the 
LMC halo, we carry out the test on each of the MACHO 
fields separately. We then evaluate and report the value 
averaged on the 30 MACHO fields. The associated disper- 
sion gives the error bars drawn in Figure [3] 

In Table [21 and Figure [3] we present the results of 
this analysis. We report the value of the parameter a 
and the KS coefficient in the three cases considered, 
^ S MW' ^ S LMC ano - ks a respectively, as a function of the 
MACHO mass. We report the results only in the range 
0.01-1 Mq where the preferred solution are found (see be- 
low). 

When we consider the MW and the LMC halos, sep- 
arately, the solutions with the highest level of confidence 
fall in correspondence of a MACHO mass « 0.5 M Q and 
« 0.15 Mq respectively. This result is confirmed by the mi- 
crolensing rate normalised distributions in Fig. [21 where 
the profile corresponding to the case of 0.5 Mq MW lenses 
is almost coincident with that of 0.2 M© LMC lenses. 

For the case of the combined populations, we find that 
the coefficient ks a presents two minima, near to the two 
values of mass found in the previous case. The absolute 
minimum, that we note to be lower than the values ob- 
tained in the case of the separate test analysis, is found at 
m = 0.15 M Q with a ~ 0.1, and gives us the parameters 
with the highest confidence level. Moreover we observe 
that the a parameter grows monotonically as a function 
of the mass from up to 1 (Fig. [SJ. 

We conclude that the statistical analysis made on the 
duration of the events gives a first suggestion that a signif- 
icant fraction of the observed events could belong to the 
LMC MACHO halo population. 

Next we consider th e issue of the spatial distribution 
of the observed events The optical depth 

profiles clearly show that the LMC halo events are char- 
acterized, with respect to LMC self-lensing and to MW 
halo ones, by a strong asymmetry with respect to the line 
of nodes (Paper II). Looking at the expected number of 
events per field, this asymmetry is somewhat weakened, 
but still present, because of the different source density as 
a function of the position. In order to get an insight in the 
more complex two dimensional situation we have to deal 
with, in Figure 0] we show the normalized number of the 
expected events (iVev), for the different lens populations 
we consider, evaluated along the axis orthogonal to the 
line of nodes passing through the LMC centre (the £ axis 
in the plot). N^v has been calculated in the LMC centre 
and in six further positions, specularly symmetric two by 
two with respect to the centre. This plot shows clearly 
that the distribution in the case of self-lensing events is 
symmetric, and moreover that outside the bar region it de- 
clines sharply. The profile for MW events presents a slight 
asymmetry with respect to the line of nodes, whereas that 
corresponding to LMC MACHO lenses have a pronounced 
asymmetric distribution. 

In Paper II we have addressed the question whether 
the observed asymmetrical distribution of the detected 
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Fig. 4. Normalized expected number profiles along the 
axis, pointing south-west, orthogonal to the line of nodes 
and passing through the LMC centre, for different lens 
populations: self lensing, dot dashed line; MW MACHO, 
solid line; LMC MACHO, dashed line. Values on the £ 
axis are in kpc. 
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Fig. 5. Location of the 30 MACHO fields in a reference 
frame centered in the LMC centre with the x-axis antipar- 
allel to the right ascension, the y-axis parallel to the dec- 
lination, and the z-axis (not shown) pointing towards the 
observer. The location of the 13 microlensing candidates, 
subset of the original set B of MACHO candidates we 
use in the present analysis ( Section I4.1|l . is shown. Also 
shown, the position of the line of nodes and the central 
band around it, which we exclude in the asymmetry anal- 
ysis. 



events, that goes indeed in the same sense predicted by a 
halo LMC population, does reflect the observational strat- 
egy, mainly to argue against the self lensing origin of the 
events. Here we take advantage of the knowledge of the ex- 
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Fig. 6. Probability isocontours with 34%, 68% and 90% regions for the LMC and MW dark matter halo fraction for 
four values of the MACHO mass. 



pected number of events for all the populations of interest 
to further study this issue. 

We take the line of nodes as the symmetry axis, and as 
in Paper II we bin the observed space in two regions, the 
north-east "near" and the south-west "far" one (Fig. EJ- 
We delimit them by two straight lines parallel to the line of 
nodes, each at a distance from the lat ter of sa 0.35 kpc, the 
1 a er ror in position as calculated bv lvan der Marel et alJ 
(2002). Note that this way the innermost LMC regions 
are excluded from the asymmetry analysis. Our purpose 
is to test the asymmetry with respect to the prediction 
of the different lens populations, assigning the probabil- 
ity pi that a microlensing event would fall in the first or 
second region to be proportional to the expected num- 
ber of events of the given population. As in Paper II we 
make use of the non-parametric Pearson test, that allows 
us to evaluate the probability to accept the null hypothe- 
sis that the pi match the measured distribution, for which 
we get 1 and 7 observed events in the near and far region, 
respectively. The analysis is carried out normalizing the 
number of expected events of each population to that of 
the observed events. 

The result of this analysis gives us a probability of 46% 
and 26% for a MACHO LMC and MW halo lens popula- 
tion, respectively. This outcome puts clearly in evidence 
the lack of a predominant lens population. Indeed, the ex- 



pected smoother spatial distribution of MW halo lenses is 
only poorly coherent with the observed distribution, thus 
challenging the explanation that attributes most of the 
events to this population. Rather, we find a much better 
agreement with the expected asymmetric distribution of 
LMC halo events. Finally, we remark that, contrary to 
the previous analysis based on the timescale distribution, 
this result turns out to be independent from the value of 
the MACHO mass. As a byproduct of the present analysis 
(confirming that carried out in Paper II) , we note that the 
probability to accept the hypothesis of a self-lensing origin 
on the basis of the observed spatial distribution turns out 
to be only of 19%. 

Both these analysis, carried out working on normal- 
ized distributions, i.e. independent from the actual halo 
fraction, indicate that a large fraction of the lenses could 
indeed belong to the LMC dark matter halo. 

4.3. MW and LMC: two different halo fractions? 

The previous analysis provided us with two important 
clues both showing that a significant fraction of the events 
detected by the MACHO collaboration could be part of 
the LMC dark matter halo. In order to reconcile this re- 
sult with the predicted number of events (Table we 
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now drop the hypothesis of equal halo MACHO fractions, 
in the MW and the LMC halo. 

We start by evaluating the likelihood function 



L (/mw> /lmc) 

exp(-A cxp ) Yl 



E8(T Ei) §-(T Ei) 



(J) 



where /]y[W, /LMC are * ne na ^° fractions for the MW 
and the LMC respectively. For both iV eX p, the expected 
number of events, and the differential rate gj^, we sum 
over all the lens populations (including the stellar ones), 
multiplying the MACHO contributions for the appropri- 
ate halo fraction. The product runs over the A D b s ob- 
served events. By Bayesian inversion, we use a flat prior 
probability, it is then possible to compute the probabil- 
ity distribution for the halo fractions given the observed 
events, P (/j^w, /lmc) - Note that we are now taking 
the MACHO mass as a parameter, equal for both halo 
populations. 

In Figure HO we show, for four values of the MACHO 
mass, the 2-dimensional probability isocontour for the two 
halo fractions. Eventually, after marginalisation over one 
variable with respect to the other, we get the results for 
the two halo fractions as a function of the MACHO mass 
shown in Figure 
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Fig. 7. MW and LMC dark matter halo fraction, median 
value with 68% CL errors, as a function of the MACHO 
mass. 

The main outcome of this analysis is that in a sig- 
nificant range of the MACHO mass, the LMC and MW 
dark matter halo fractions are not expected to be equal. 
In particular, from 0.1 up to 0.3 Mq, both a large value 
for /lmc an d a small one for /mw are expected. 

This behaviour is mainly due to the observed 
timescales. It is therefore not surprising that this re- 



sult is coherent with that obtained with the KS test, 
where we have found a preferred value of the mass of 
about 0.2 Mq, with a significant expected contribution 
from LMC MACHO halo lenses. The likelihood analysis 
gives, for m = 0.2 Mq, / mw ~ 4% and /lmC ~ 80 %- 
At face value, given the number of expected events, we 
get to the overall prediction of about 6-7 events to be at- 
tributed to the LMC dark matter halo lens population, 
2-3 to the MW halo one, to be looked for among those 
with the shorter durations, still allowing for a contribu- 
tion, about 2-3 events, distributed among self-lcnsing and 
LMC stellar halo lenses. 

For larger values of the mass, the LMC MACHO halo 
fraction turns out to be almost degenerate though com- 
patible with zero, and, for 0.5 M , we recove r the result of 
a MW MACHO halo fraction of about 20% l|Alcock et alJ 
l2000b|) . This is again coherent with the issue of the KS 
test. 

We stress that the outcome of this analysis has to be 
looked at together with the outcomes of the previous anal- 
ysis on the timescale and the spatial distributions of the 
observed events. Overall, they indicate that i) a sizeable 
fraction of the observed events shows characteristics in ac- 
cord with those expected for a MACHO LMC halo popula- 
tion; ii) such a contribution may be expected by dropping 
the hypothesis that the halo fractions in form of MACHOs 
in the MW and the LMC dark matter halo are equal. 

As a last point we discuss the dependence of our results 
from the choice of the LMC dark matter halo parameters, 
the central density and the truncation radius (Sect. |2J). 
We find that any variation of the parameters in a reason- 
able range around their fiducial values do not change, at 
least qualitatively, our results. As for the central density, 
any variation downward (upward) is linearly related with 
a corresponding change in the total number of expected 
events. This implies a scaling respectively upward (down- 
ward) for the halo fraction without affecting the main con- 
clusion on the contribution of the LMC halo. As for the 
truncation radius (Rt), in first approximation the situa- 
tion is similar, as, roughly, a smaller (larger) value of Rt 
give rises to a smaller (larger) total LMC mass. However, 
the issue is slightly more subtle. First, the problem is not 
symmetric with respect to the choice of our fiducial value, 
Rt — 15 kpc. Second, different choices for R t enhance 
different spatial distribution for the number of expected 
events. In particular it results that: i) the decrease in the 
number of expected events, for values of Rt smaller than 
the fiducial value, is relatively larger than the increase ob- 
tained by choosing larger values; ii) the spatial asymme- 
try of MACHO LMC halo events is enhanced for smaller 
value of Rt- (Both effects are easily explained as most of 
the lenses are expected to be located in the innermost 
LMC region.) We have tested our results with 4 values of 
Rt, respectively smaller (larger) of our fiducial value by 
20% and 40%. The relative average decrease in the num- 
ber of expected events for R t = 9 and 12 kpc with respect 
to the Rt = 15 kpc case is ~ 20% and 8.7%, whereas the 
relative increase for R t = 18 and 21 kpc is ~ 7.4% and 
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14% respectively. As for the spatial distribution, we char- 
acterise the asymmetry by the relative difference between 
the number of expected events evaluated at the extremes 
of the £ axis (as defined in Fig. |3J| £ = 1 kpc and £ = — 1 
kpc. With respect to the fiducial case, where this turns 
out to be of 47%, for R t = 9 and 12 kpc we find a relative 
increase of 11% and 4%, whereas for R t = 18 and 21 kpc 
a relative decrease of 2% and 6%, respectively. Overall, 
smaller values of Rt strengthen our conclusions. 

5. Conclusions 

In this paper we have addressed the issue of the inter- 
pretation of the microlensing results toward the LMC. In 
particular, starting from the sample of microlensing can- 
didates reported by the MACHO collaboration, we have 
discussed the contingent contribution of a lens component 
belonging to the LMC dark matter halo besides that of 
the MW one. As a main result of the present analysis we 
show that a sizeable fraction of the observed events, up to 
about half of the total, could indeed be part of the first 
component. 

We summarise our analysis as follows. First, we have 
compared the observed timescales with those expected for 
the two different MACHO populations, the MW and the 
LMC one. As a result we have shown that the preferred 
values for the MACHO mass are about 0.5 and 0.2 M Q re- 
spectively and, through a KS test, that the latter solution 
is preferred. Second, we have studied the spatial distribu- 
tion of the observed events, recalling that, because of the 
inclination of the LMC disc with respect to our line of 
sight, an asymmetry is expected for LMC halo events. As 
a result we have shown that, independently from the value 
of the MACHO mass, the observed distribution matches 
better that expected for a LMC halo population with re- 
spect to that of a MW halo population. Overall, these 
are clues suggesting the presence of a significant MACHO 
LMC halo population among the observed events. 

The extremely larger value of the overall MW halo 
mass with respect to the LMC one implies that generally 
one can safely ignore the LMC halo component. In order to 
explain such a large contribution of the latter, one way out 
is to consider that the halo fractions in form of MACHOs 
of the two components, the MW and the LMC halos, may 
be different. Coherently with the timescale analysis, this 
issue turns out to be strongly dependent on the value of 
the MACHO mass. 

In order to get to more quantitative results we have 
tested this hypothesis through a likelihood analysis. 
Eventually we have shown that for a large range of mass 
values a different (and larger) fraction for the LMC halo 
with respect to the MW one is indeed expected. In par- 
ticular, for MACHO mass of ~ 0.2 M Q we evaluate a high 
halo fraction for the LMC, ~ 80%, together with a small 
one for the MW, 5%, thus implying that about half of 
the observed events should belong to the LMC dark mat- 
ter halo. On the other hand, for MACHOs of ~ 0.5 M Q 
we recover the well known result of a MW halo fraction 



~ 20% with a (possibly) negligible contribution from the 
LMC dark matter halo. 

A possible explanation to the origin of different halo 
fractions could come from the different formation history 
of the two galaxies, or, more simply, could be related to the 
fact that one observes all the LMC halo but, practically, 
only a line of sight through the Galactic halo. 

These conclusions should be taken cum grano salis. 
The overall implicit assumption is the validity of the 
MACHO results, whereas they are actually challenged by 
the EROS collaboration. An intrinsic limit is then due to 
the sm all statistic at dispo sal. The SuperMACHO collabo- 
ration l|Becker et alJl200^l is expected to provide a larger 
sample of candidates spread over a much larger field of 
view and this should allow to put firmer constraints on 
this problem. Furthermore, the model issue, in particular 
for the LMC components, is still a matter of debate. Our 
analysis shows, however, that it is in principle possible to 
characterize and distinguish the two halo lens populations 
and, moreover, challenge the usual implicit assumption of 
equal halo fractions in form of MACHO for both the MW 
the LMC dark matter halos. 
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